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Abstract. We present measurements of the transverse momentum spectra, yield and 
(pt) systematics for Kg, A and A in p + p collisions at y/s = 200 GeV. We show a 
dependence of the (px) with event multiplicity and infer that this is consistent with a 
mini-jet dominated particle production mechanism. These observations are compared 
to available data from p + p experiments as well as to pQCD theoretical predictions. 



1. Introduction 

Particles which contain strange quarks are valuable probes of the dynamics of p + p 
collisions as constituent strange quarks are not present in the initial colliding nuclei. 
The enhancement of the strange particle yield from p+p collisions to heavy ion collisions 
has been suggested as a possible Quark Gluon Plasma signature [T] . 
We report preliminary results from the 2001/2002 p+p run by the STAR experiment. 
The main focus of this paper is on presenting a high statistics measurement of Kg, A 
and A at yfs = 200 GeV and obtaining the yield and (pt) for each species. A discussion 
of appropriate parameterizations to the particle spectra will ensue. We will show that 
the widely used power- law extrapolation in p + p collisions is not sufficient to obtain the 
best x 2 results for fits to the strange particle spectra and we will consider alternatives [2]. 
Furthermore, we will study the dependency of (px) and spectral shape as a function of 
charged particle event multiplicity(A c /j), an effect that has been observed previously in 
high energy p + p collisions. Several authors ascribe this phenomenon to pQCD mini- 
jets fragmenting into hadronic final states [HUE]. However, we show that present string 
models incorporating 'semi-hard' pQCD processes, such as PYTHIA, do not describe 
the STAR data well without significant parameter tuning 5j. 
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2. Analysis 

2.1. Event selection and Pile-up 

The data were reconstructed using the STAR detector system which is described in 
more detail elsewhere [01 Ej . The main tracking detector used in this analysis is the 
Time Projection Chamber (TPC) covering the full acceptance in azimuth and a large 
pseudo-rapidity coverage (| n |< 1.5). 

A total of 14 million non-singly diffractive (NSD) events were triggered with the STAR 
beam-beam counters (BBC) requiring two coincident charged tracks at forward rapidity 
(3.5 <| 7] |< 5.0). Due to the particulary low track multiplicity environment in p+p 
collisions a special vertex finder was developed. A Monte-Carlo study showed that only 
76% of primary vertices are found correctly; from the remainder, 14% are lost and 10% 
are badly reconstructed. A total of 11 million events passed the event level selection 
criteria requiring a valid primary vertex within 100cm along the beam-line from the 
center of the TPC detector. An additional multiplicity dependant vertex efficiency 
factor e vtx (N ch ) must therefore be applied to the particle spectra. 

The high luminosity of the RHIC proton beams has the unfortunate effect of producing 
several collisions during the read-out time of the TPC, a condition known as "pile-up". 
These pile-up events cannot be identified separately, however the triggered vertex can be 
identified by allowing only tracks that match to a very fast scintillating Central Trigger 
Barrel(CTB) detector to be used in the vertex reconstruction. 

2.2. Particle reconstruction 

The strange particles were identified from their weak decay to charged daughter particles 
via topological identification algorithm. The following decay channels were analyzed: 
A -> p + vr-(b.r. 63.9%) ,A -> p + vr+(b.r. 63.9%), K° s -> tt+ + tt - (b.r. 68.6%). 
Particle identification of daughters was achieved by requiring the dE/dx to fall within 
the 3cx-bands of the theoretical Bethe-Bloch parameterizations. Further background in 
the invariant mass was removed by applying topological cuts to the decay geometry, 
i.e. requiring the parent to originate from within 2cm of the primary vertex and have a 
minimum decay length of 2cm. These cuts also ensure that the particle decay is matched 
to the correct vertex with a high probability. 

Corrections for particle reconstruction efficiency were obtained through a Monte-Carlo 
based method of embedding simulated particle decays into real events and comparing the 
number of simulated and reconstructed particles in each p^-bin. The A and A spectra 
were corrected for feed-down from H (H~ — > A + 7r~). Finally, the corrected transverse 
momentum spectra have to be multiplied by the multiplicity dependant primary vertex 
efficiency e vtx (N c h) to account for yield losses. 
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Figure 1. Minimum-bias, non-feeddown corrected spectra for Kg (left) and A (right) 
with (| y |< 0.5) from p + p at ^/s = 200GeV. The errors in the plot are statistical 
only. The fit is composite and is described in the text. 



3. Transverse Momentum Spectra 

Non-feeddown corrected spectra for Kg, A and A are shown in Figure H The particle 
acceptance at mid-rapidity (| y |< 0.5) in the TPC starts at transverse momentum 0.2 
GeV/c for Kg and 0.3 GeV/c for A and A. The spectra are parameterized to allow 
extrapolation to low pt and to extract (px) and yield (dN/dy). In contrast to previous 
p + p experiments j2], which used either a single exponential function in transverse mass 
or power-law functions, we found that a combination of these functions is more effective 
in fitting the STAR data. The x 2 results from different fit-functions are presented in 
Table [T] Composite fits, represented in equation 1 (for Kg) and equation 2 (for A and 
A) yield the lowest x 2 an d were used to extract the final values for (px) and yield as 
shown in Table El 
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Table 1. A summary of \ 2 values for different fit-functions to the particle spectra 
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The systematic errors are dominated by the different possible fit parameterizations 
and amount to 10% for (px) and 15% on the yield. The remaining systematical 
uncertainties are from the background and efficiency corrections. Table El compares 
the STAR results to the measurements by UA5 [Hj. The UA5 measurements were made 
over a larger rapidity interval (| y |< 2) and have been scaled to match the STAR 
acceptance by the use of a PYTHIA simulation. The A/A ratio is 0.88±0.09 and is flat 
as a function of px, consistent with a nearly net-baryon free environment at mid-rapidity. 
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Table 2. A summary of yield and (px) for measured strange particles. The errors 
quoted are statistical only. 
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Table 3. A summary of yield and (px) for and A, A (feed down corrected) and Kg 
measured by the UA5 [EJ ■ UA5 yields have been scaled to the STAR measured rapidity 
interval using a Pythia simulation 



4. Comparison to PYTHIA 

Figure El shows the comparison of the STAR spectra to the spectra from the default 
setting of PYTHIA event generator (version 6.22, MSEL1) which includes hard processes 
and a modelling of soft processes, through a pQCD model. 

Even at px > 1.5, a regime were pQCD should give a better description, the 
differences between the model and the data are considerably large. Further studies 
are planned in order to optimally tune the model to match the STAR strangeness data. 

5. (px) vs event multiplicity 

The minimum bias event sample was split into 6 event classes according to mean charged 
particle multiplicity per unit rj ((dN c h/drj)). A similar analysis has been performed in 
previous p + p experiments at different collision energies jSHH]- High multiplicity events 
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are particularly interesting as they are likely dominated by contributions from large 
momentum transfer parton-parton collisions. 
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Figure 3. (px) vs. N c h for Kg (left) and A (right) compared to PYTHIA (triangles) 
(ver 6.22, MSEL1). The errors on data are both statistical and systematic. 



Figure El presents a comparison of the (px) vs (dN c h/drj) for A and Kg with the 
corresponding PYTHIA calculations. A rise in (px) with increasing N c h is observed and 
is stronger for the heavier particle (A). Several authors have attributed these phenomena 
to the increased number of fragmenting mini-jets in the high- multiplicity events [TTlll2j . 
In fact, the HIJING model even makes a prediction for the expected number of mini-jets 
contributing to each multiplicity [TO] . 
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It is also interesting that PYTHIA under-predicts the magnitude and correlation 
strength of (px) with multiplicity. Similar comparisons for inclusive charged particles 
(Pt) have been recently published by the CDF Collaboration and demonstrate similar 
discrepancies |T3] . 

6. (px) vs particle mass 

In Figure 0] the (px) is shown for all particle species measured in p + p collisions in 
STAR as a function of particle mass. The black line is an empirical parameterization 
to ISR data for vr,K,p only, at y/s = 25 GeV. There is very little energy dependance for 
the lower mass particles between the ISR and RHIC energies, whereas the higher mass 
particles show a clear breaking from the parameterization. The dependance of (px) on 
particle mass has been explained by mini-jet production in p + p and p + p collisions 
and a stronger contribution for higher mass particles is expected jlj. 




Figure 4. (px) vs particle mass, where the black curve represents the ISR 
parametrization from tt, K ,p at y/s = 25 GeV. 



7. Summary 

The STAR experiment has made the first high statistics measurement of mid-rapidity 
Kg, A and A in p + p collisions at yfs = 200 GeV. The results agree with those made by 
the UA5 collaboration for p+p at yfs = 200 GeV. The ratio of A/A suggests a small net 
baryon number at mid-rapidity. We have show that the shape of the particle spectra 
is best described with a two-component fit function accounting for both the soft and 
hard parts of the spectra. Furthermore, we have undertaken studies to understand the 
change in shape and (px) of the spectra with increased event multiplicity. We compared 
our results to pQCD inspired models which infer mini-jets and string fragmentation. We 
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have shown that the current default version of PYTHIA is not well tuned for strange 
particles at this collision energy. 
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2. Analysis 

2.1. Event selection and Pile-up 

The data were reconstructed using the STAR detector system which is described in 
more detail elsewhere [6, 7]. The main tracking detector used in this analysis is the 
Time Projection Chamber (TPC) covering the full acceptance in azimuth and a large 
pseudo-rapidity coverage (| r\ |< 1.5). 

A total of 14 million non-singly diffractive (NSD) events were triggered with the STAR 
beam-beam counters (BBC) requiring two coincident charged tracks at forward rapidity 
(3.5 <| i] |< 5.0). Due to the particulary low track multiplicity environment in p+p 
collisions a special vertex finder was developed. A Monte-Carlo study showed that only 
76% of primary vertices are found correctly; from the remainder, 14% are lost and 10% 
are badly reconstructed. A total of 11 million events passed the event level selection 
criteria requiring a valid primary vertex within 100cm along the beam-line from the 
center of the TPC detector. An additional multiplicity dependant vertex efficiency 
factor e v tx{N c h) must therefore be applied to the particle spectra. 

The high luminosity of the RHIC proton beams has the unfortunate effect of producing 
several collisions during the read-out time of the TPC, a condition known as "pile-up". 
These pile-up events cannot be identified separately, however the triggered vertex can be 
identified by allowing only tracks that match to a very fast scintillating Central Trigger 
Barrel(CTB) detector to be used in the vertex reconstruction. 

2.2. Particle reconstruction 

The strange particles were identified from their weak decay to charged daughter particles 
via topological identification algorithm. The following decay channels were analyzed: 
A -> p + 7r-(b.r. 63.9%) ,A -> p + 7r+(b.r. 63.9%), K° -> tt+ + tt" (b.r. 68.6%). 
Particle identification of daughters was achieved by requiring the dE/dx to fall within 
the 3<r-bands of the theoretical Bethe-Bloch parameterizations. Further background in 
the invariant mass was removed by applying topological cuts to the decay geometry, 
i.e. requiring the parent to originate from within 2cm of the primary vertex and have a 
minimum decay length of 2cm. These cuts also ensure that the particle decay is matched 
to the correct vertex with a high probability. 

Corrections for particle reconstruction efficiency were obtained through a Monte-Carlo 
based method of embedding simulated particle decays into real events and comparing the 
number of simulated and reconstructed particles in each p^-bin. The A and A spectra 
were corrected for feed-down from S (S _ — > A + tt~). Finally, the corrected transverse 
momentum spectra have to be multiplied by the multiplicity dependant primary vertex 
efficiency e vtx (N ch ) to account for yield losses. 
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Figure 1. Minimum-bias, non-feeddown corrected spectra for Kg (left) and A (right) 
with (| y |< 0.5) from p + p at y/s = 200GeV. The errors in the plot are statistical 
only. The fit is composite and is described in the text. 



3. Transverse Momentum Spectra 

Non-feeddown corrected spectra for Kg, A and A are shown in Figure 1. The particle 
acceptance at mid-rapidity (| y |< 0.5) in the TPC starts at transverse momentum 0.2 
GeV/c for Kg and 0.3 GeV/c for A and A. The spectra are parameterized to allow 
extrapolation to low pt and to extract (px) and yield (dN/dy). In contrast to previous 
p + p experiments [2], which used either a single exponential function in transverse mass 
or power-law functions, we found that a combination of these functions is more effective 
in fitting the STAR data. The x 2 results from different fit-functions are presented in 
Table 1. Composite fits, represented in equation 1 (for Kg) and equation 2 (for A and 
A) yield the lowest x 2 an d were used to extract the final values for (px) and yield as 
shown in Table 2. 
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Table 1. A summary of \ 2 values for different fit-functions to the particle spectra 
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The systematic errors are dominated by the different possible fit parameterizations 
and amount to 10% for (px) and 15% on the yield. The remaining systematical 
uncertainties are from the background and efficiency corrections. Table 3 compares 
the STAR results to the measurements by UA5 [8]. The UA5 measurements were made 
over a larger rapidity interval (| y |< 2) and have been scaled to match the STAR 
acceptance by the use of a PYTHIA simulation. The A/A ratio is 0.88±0.09 and is flat 
as a function of px, consistent with a nearly net-baryon free environment at mid-rapidity. 
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Table 2. A summary of yield and (pt) for measured strange particles. The errors 
quoted are statistical only. 
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Table 3. A summary of yield and (px) for and A, A (feed down corrected) and Kg 
measured by the UA5 [8] . UA5 yields have been scaled to the STAR measured rapidity 
interval using a Pythia simulation 



4. Comparison to PYTHIA 

Figure 2 shows the comparison of the STAR spectra to the spectra from the default 
setting of PYTHIA event generator (version 6.22, MSEL1) which includes hard processes 
and a modelling of soft processes, through a pQCD model. 

Even at pt > 1.5, a regime were pQCD should give a better description, the 
differences between the model and the data are considerably large. Further studies 
are planned in order to optimally tune the model to match the STAR strangeness data. 

5. (px) vs event multiplicity 

The minimum bias event sample was split into 6 event classes according to mean charged 
particle multiplicity per unit rj ((dN c h/ 'drj)). A similar analysis has been performed in 
previous p + p experiments at different collision energies [2, 9]. High multiplicity events 
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Figure 2. Kg (left) and A (right) particle spectra (circles) compared to PYTHIA 
(line) (ver 6.22, MSEL1) 



are particularly interesting as they are likely dominated by contributions from large 
momentum transfer parton-parton collisions. 
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Figure 3. (p T ) vs. N ch for Kg 1 (left) and A (right) compared to PYTHIA (triangles) 
(ver 6.22, MSEL1). The errors on data are both statistical and systematic. 



Figure 3 presents a comparison of the (px) vs (dN c h/di]) for A and Kg with the 
corresponding PYTHIA calculations. A rise in (px) with increasing N c h is observed and 
is stronger for the heavier particle (A). Several authors have attributed these phenomena 
to the increased number of fragmenting mini-jets in the high-multiplicity events [11, 12]. 
In fact, the HIJING model even makes a prediction for the expected number of mini-jets 
contributing to each multiplicity [10]. 
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It is also interesting that PYTHIA under-predicts the magnitude and correlation 
strength of (px) with multiplicity. Similar comparisons for inclusive charged particles 
(px) have been recently published by the CDF Collaboration and demonstrate similar 
discrepancies [13]. 

6. (p T ) vs particle mass 

In Figure 4 the (px) is shown for all particle species measured in p + p collisions in 
STAR as a function of particle mass. The black line is an empirical parameterization 
to ISR data for 7r,K,p only, at ^/s = 25 GeV. There is very little energy dependance for 
the lower mass particles between the ISR and RHIC energies, whereas the higher mass 
particles show a clear breaking from the parameterization. The dependance of (pt) on 
particle mass has been explained by mini-jet production in p + p and p + p collisions 
and a stronger contribution for higher mass particles is expected [4]. 



aj 4 
v i.t 

1.2 

1 

0.8 
0.6 
0.4 
0.2 

°0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

Mass (GeV/c 2 ) 

Figure 4. (px) vs particle mass, where the black curve represents the ISR 
parametrization from 7r, K ,p at yfs = 25 GeV. 



7. Summary 

The STAR experiment has made the first high statistics measurement of mid-rapidity 
Kg, A and A in p + p collisions at ^/s = 200 GeV. The results agree with those made by 
the UA5 collaboration for p + p at ^/s = 200 GeV. The ratio of A/ A suggests a small net 
baryon number at mid-rapidity. We have show that the shape of the particle spectra 
is best described with a two-component fit function accounting for both the soft and 
hard parts of the spectra. Furthermore, we have undertaken studies to understand the 
change in shape and (px) of the spectra with increased event multiplicity. We compared 
our results to pQCD inspired models which infer mini-jets and string fragmentation. We 
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have shown that the current default version of PYTHIA is not well tuned for strange 
particles at this collision energy. 
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